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We present a model for structure formation, melting, and optical properties of gold/DNA 
nanocomposites. These composites consist of a collection of gold nanoparticles (of radius 50 nm 
or less) which are bound together by links made up of DNA strands. In our structural model, the 
nanocomposite forms from a series of Monte Carlo steps, each involving reaction-limited cluster- 
cluster aggregation (RLCA) followed by dehybridization of the DNA links. These links form with a 
probability p e // which depends on temperature and particle radius a. The final structure depends 
on the number of monomers (i. e. gold nanoparticles) N m , T, and the relaxation time. At low 
temperature, the model results in an RLCA cluster. But after a long enough relaxation time, the 
nanocomposite reduces to a compact, non-fractal cluster. We calculate the optical properties of 
the resulting aggregates using the Discrete Dipole Approximation. Despite the restructuring, the 
melting transition (as seen in the extinction coefficient at wavelength 520 nm) remains sharp, and 
the melting temperature Tm increases with increasing a as found in our previous percolation model. 
However, restructuring increases the corresponding link fraction at melting to a value well above the 
percolation threshold. Our calculated extinction cross section agrees qualitatively with experiments 
on gold/DNA composites. It also shows a characteristic "rebound effect," resulting from incomplete 
relaxation, which has also been seen in some experiments. We discuss briefly how our results relate 
to a possible sol-gel transition in these aggregates. 



I. INTRODUCTION 

The optical properties of metallic nanoparticles have 
been investigated intensively for many years [1,2]. Re- 
cently, this work has extended to so-called functional 
metallic nanoparticles, which may have a variety of novel 
and useful optical and mechanical properties [3] . Among 
these, there has been particular interest in the DNA- 
modified gold nanoparticle (gold/DNA nanocomposite) 
system. This is a material consisting of gold nanopar- 
ticles to which specific kinds of organic molecules (e. 
g., noncomplementary oligonucleotides capped with thiol 
groups) can be attached. These materials can be pro- 
duced in a variety of structures using a strategy of 
nanoparticle self-assembly [4-8] . They may also be useful 
for selective biological detection, by making use of the op- 
tical and electrical sensitivity of their aggregates [9-14] . 

Numerical model calculations of the optical proper- 
ties of DNA modified gold nanoparticle aggregates show 
general agreement with experiments. For example, both 
experiments and calculations show that (i) for isolated 
gold nanoparticles in suspension, there is a strong sur- 
face plasmon absorption in the visible, arising from oscil- 
lations of electronic charge in the gold nanoparticles; and 
(ii) this absorption maximum broadens and red-shifts 
when the cluster radius becomes comparable to the wave- 
length [15-17]. 

Although the DNA molecules absorb primarily in the 
ultraviolet, they nonetheless play a central role in the op- 



tical properties of these aggregates in the visible, because 
they strongly influence the structure of the gold/DNA ag- 
gregates. For example, the DNA tends to form multiple 
links between individual gold nanoparticles. These mul- 
tiple links appear to account for some key structural fea- 
tures of the aggregates [18-20]. In particular, the melting 
transition for a gold nanoparticle aggregate has a much 
narrower temperature width, and occurs at a substan- 
tially higher temperature T, than that of a single DNA 
duplex [9,21]. In addition, the presence of multiple links 
leads to a natural explanation for the dependence of the 
aggregate melting temperature, denoted Tm on particle 
size [20,22]. 

In the present work, we extend our previous model 
calculations [18,19] to take into account the dynamics of 
aggregate formation, and how these dynamics affect the 
aggregate optical properties. In our previous work, the 
melting of the aggregates was treated using a purely sta- 
tistical criterion. Specifically, the aggregates were formed 
by removing DNA links between gold nanoparticles on 
a simple cubic lattice, with a suitable, temperature- 
dependent probability. In the limit of a large aggregate, 
the melting transition occurs, in this model, when the 
fraction of links falls below the percolation threshold p c 
for the lattice considered. The calculated aggregate opti- 
cal properties are found to change dramatically when this 
threshold is passed, in good agreement with experiment. 

To improve on this approach, we describe below a 
model for structure formation which starts from isolated 
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gold nanoparticles. Our model takes into account two im- 
portant features of the structure formation: the diffusion 
of nanoparticles, or clusters of nanoparticles, through the 
solvent to form a cluster, and the chemical reaction be- 
tween DNA chains which produces the links between the 
nanoparticles. 

Our structural model leads to a wide range of possi- 
ble aggregate morphologies, depending on the temper- 
ature. Corresponding to these morphologies is a broad 
range of possible optical properties. In this paper, we 
will present numerical results for both the structural and 
optical properties of these nanocomposites over a typical 
range of parameters. For physically reasonable param- 
eters, our numerical results are in good agreement with 
experiment. We will present a qualitative interpretation 
of these results, and compare them to our earlier, purely 
percolation model. 

The remainder of this paper is organized as follows. 
In Section II, we describe our structural model for 
gold/DNA nanocomposites, and explain how it is imple- 
mented numerically. In Section III, we review the Dis- 
crete Dipole Approximation, which is the method used 
to calculate the optical properties of this system; we also 
discuss various technical details needed to treat the ir- 
regular clusters which emerge from the structural model 
of Section II. In Section IV, we present our numerical re- 
sults for both the structural and the optical properties of 
the aggregates. Finally, in Section V, we summarize our 
results, interpret them in terms of the expected behavior 
of typical gold/DNA nanocomposites, compare our re- 
sults to available experiments, and discuss their possible 
implications for future work. 

II. STRUCTURAL MODEL AND ITS 
NUMERICAL IMPLEMENTATION 

In this Section, we describe our structural model for 
the formation of gold/DNA aggregates. We start by de- 
scribing the expected aggregation behavior at low T. Fol- 
lowing this, we present our full structural model and its 
numerical implementation for arbitrary T. 

A. Aggregation at Low T. 

Before discussing the aggregation of gold/DNA 
nanocomposites at low T, we first consider the aggrega- 
tion behavior of other typical colloids, i. e. suspensions 
of small solid particles in a liquid solvent. As the individ- 
ual colloidal particles stick together to become clusters, 
these clusters themselves diffuse through the solvent, and 
continue to collide and aggregate. This behavior is an ex- 
ample of cluster-cluster aggregation. 

If bonding between two colloidal particles is irre- 
versible, the final aggregate generally results from one 



of two processes: reaction-limited aggregation (RLA) or 
diffusion- limited aggregation(DLA). If there is a repulsive 
energy barrier between two approaching colloidal parti- 
cles, the resulting process is expected to be a reaction- 
limited aggregation, because the reaction barrier is the 
limiting step in cluster growth [23]. By contrast, if there 
were no barrier, the aggregation would be dominated by 
diffusion processes, and the resulting clusters should ex- 
hibit the features of DLA. DLA clusters, like RLA clus- 
ters, are fractals, but they have a substantially lower df 
than the RLA clusters. For forming both types of fractal 
clusters, however, the irreversibility of binding is impor- 
tant. Without irreversibility, the final cluster is likely to 
become compact and non- fractal [24,25]. 

We now provide a possible justification for considering 
the growth of the colloidal cluster in a dilute solution of 
DNA-modified gold particles at very low T as reaction- 
limited cluster-cluster aggregation (RLCA). Here, we 
should note that for the DNA-modified gold nanoparticle 
system, the mechanism for binding two DNA-modified 
gold nanoparticles differs from the ordinary colloid ag- 
gregation process discussed above, because binding can 
occur only if the DNA hybridizes, i. e., if the two DNA 
single strands on different gold particles and one linker 
single DNA strand undergo a reversible chemical reac- 
tion to form a double strand. However, for two reasons, 
RLCA may still be a plausible growth mechanism for the 
gold cluster aggregation in a dilute solution at very low 
T. First, at the very low T, once the DNA hybridization 
occurs and the two gold particles do stick together, they 
rarely unbind, since thermal fluctuation cannot provide 
enough energy to break them apart. Thus, the aggre- 
gation precess is irreversible. Secondly, in case of DNA 
hybridization, this reaction barrier can actually be exper- 
imentally observed [26]. Because of this barrier, nanopar- 
ticles must collide numerous times before two particles 
can stick together, since DNA hybridization takes a fi- 
nite time. 

As T increases, the above argument, based on the ir- 
reversible binding, is no longer valid, since DNA dehy- 
bridization may easily take place. Thus, the restructur- 
ing of clusters is crucial for understanding the cluster 
morphology and the melting transition. We now describe 
a structural model which includes the effects of this re- 
structuring on the cluster morphology. 

B. Model for Structure Formation at General 
Temperature 

1. Description of the Model 

At finite T, cluster restructuring is sensitive to the rel- 
ative magnitudes of two times scales, which we denote 
Tund and Tdehyb- Tu-nd is the time elapsed when two clus- 
ters meet and attempt to form a link, while Tdehyb is the 
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time needed for a a double DNA strand to dchybridizc 
into two single DNA strands under the influence of ther- 
mal fluctuations. At low T, Tdehyb ^> Tbmd [26]. Hence, 
there will be many cluster-cluster binding events in the 
time required for a one dehybridization. In the present 
model, we will assume that Tdehyb ^> Tu n d at all T. With 
this assumption, we can consider aggregation and dehy- 
bridization processes separately. That is, we can estab- 
lish a model in which, first the monomers undergo aggre- 
gation and form one single large cluster, and, next, the 
cluster breaks into smaller parts via a suitable cluster 
unbinding process due to DNA dehybridization [27]. 

Our aggregation procedure is carried out by a model of 
RLCA [28], as justified in the previous subsection. In this 
algorithm, the system is assumed to be made up initially 
of a large number N m of "monomers" (each consisting of 
a single gold nanosphcrc of radius a). In the first step 
of the aggregation, we choose two monomers at random, 
then place them at two different points chosen randomly 
on the sites of a simple cubic lattice of lattice constant I 
and edge L = N£ in d dimensions (d = 3 in all of our sim- 
ulations), with free boundary conditions. If the particles 
happen to be placed on adjacent sites, they are assumed 
to form a two-particle cluster, and are removed from the 
lattice. If they are placed on non-adjacent sites, the pro- 
cedure is repeated until they do form a cluster. In the 
next step, two of the N m — 1 clusters, chosen at random, 
are placed in random, but non-overlapping positions and 
random orientations on the lattice. If they are adjacent, 
the two clusters are assumed to merge and form a larger 
cluster; otherwise, the procedure is repeated until the 
number of clusters is again reduced by one. In the n th 
step, the same procedure is carried with two clusters ran- 
domly chosen from the now N m — n + 1 clusters. Even- 
tually, this procedure leads to the formation of a single 
large cluster. 

Once our aggregation procedure is finished, the next 
step in the procedure is to simulate cluster unbinding 
due to DNA dehybridization. Here we use a percolation 
algorithm, where the key parameter is p e ff, the prob- 
ability that a given link between two monomers in the 
cluster remains the same without breaking. Note that 
this p e f f is obviously related to DNA hybridization, and 
thus depends on temperature. (The exact relation of 
p e f f to DNA hybridization is discussed in the next sub- 
section.) With probability l—p e ff, we randomly remove 
links in the cluster, which has formed from the aggre- 
gation procedure. After random removal of links, the 
cluster may separate several clusters, and thus we use a 
simple computer algorithm in order to identify the sep- 
arate clusters [30] . If the resulting aggregate consists of 
two or more clusters, these aggregates are placed in ran- 
dom, but non-overlapping, positions and orientations on 
the large cubic lattice. This configuration will be used 
for the calculation of optical properties. 

The two procedures described above are defined as one 



single Monte Carlo (MC) "step." To generate the final 
configuration, we carry out a series of Nmc MC steps, 
each involving aggregation and cluster unbinding. The 
end result is one or more final clusters. 

The structure and the number of the final cluster or 
clusters depend on p e ff- If Peff = 1, the final result 
is simply a single cluster formed from the aggregation 
part of the first MC step, which for large enough N m , 
is known to be a RLCA fractal, with fractal dimension 
df = 2.1 ii d = 3 [28]. For p e ff < 1, two types of transi- 
tions can take place in this system: a "sol-gel transition" 
and a "melting transition." lip e ff is slightly less than 1, 
the restructuring of cluster does occur and tends to make 
the cluster more compact, and it may eventually become 
non- fractal at a certain p, denoted psg- We can call 
this transition as "sol-gel" transition [29]; for p > psg, 
the system will be a "gel," i. e., a fractal cluster, while 
for p < psg, it will be "sol," characterized by compact 
and non- fractal clusters. If p e ff is substantially smaller 
than 1, there are no longer large aggregation clusters. 
Thus we can define the "melting transition" as the point 
where, even in the limit of very large N m , the aggregation 
procedure described above leads only to finite clusters - 
there is no "infinite cluster." This transition occurs at a 
critical value of p e ff, denoted pu- If there is no restruc- 
turing effect, this transition is related to the bond perco- 
lation transition, and thus we can consider that pu = Pc, 
where p c is the bond percolation threshold (for example, 
p c <~ 0.25 on a 3-dimensional simple cubic lattice) [30]. 
Next, we discuss the connection between p e ff an d the 
physical parameters of the real gold/DNA nanoparticle 
system, as previously analyzed in Ref. [18]. 

2. Determination ofp e ff by DNA hybridization 

At a low temperature T, a "link" is expected to consist 
of a number, say Nd, of DNA "duplexes," i. e., of pairs of 
DNA strands connected to form a molecule. In actuality, 
there is a linker molecule which emerges from solution 
to connect two DNA single strands, each on a different 
nanoparticle. Even in this case, however, the melting 
condition can be reduced to that used here, without loss 
of accuracy, as will be shown in Ref. [19]. Also, there are 
actually two chemically distinct DNA single strands (de- 
noted A and B); each gold particle has either all A or all B 
single strands attached, and the linker molecule can con- 
nect only A and B single strands. The existence of two 
species can be disregarded by symmetry when the the 
concentrations of the two species in solution are equal; 
experiments have been carried out only under these con- 
ditions [4,9,17,20,22]. 

Thus, we simply assume that each DNA duplex con- 
sists of one double strand D, made up of two short single 
strands S (each having 12-14 DNA base pairs). To de- 
scribe DNA dehybridization, we adopt a simple two-state 
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model [26,31]. There also exists a more elaborate the- 
ory which can account for many detailed features of the 
gold/DNA nanocomposite, including the dependence of 
the melting transition temperature Tm on the salt con- 
centration of the solvent [20]. However, as mentioned in 
Introduction, both theories reached the same qualitative 
conclusion about the melting transition, namely that the 
observed sharp melting transition and the melting tem- 
perature dependence on the system size originate from 
the presence of multiple DNA links between each pair of 
nanoparticles. Moreover, the simple model we use here 
can be easily extended to the similar systems discussed 
in Rcf. [10,32-34]. 

In a simple two-state model, the relative proportion of 
D and S is determined by the chemical reaction 

S + S^D. (1) 

The chemical equilibrium condition corresponding to eq. 
(1) is 

(i-p(r)) 2 _ k(t) 

P(T) C T ' (2) 

where p(T) is the fraction of the single DNA strands 
which form double strands by the reaction (1) at tem- 
perature T, K(T) is a chemical equilibrium constant, and 
Ct is the molar concentration of single DNA strands in 
in the sample, Since < p(T) < 1, the physical solution 
to eq. (2) is 

p(T) = 1 + 1 (K' - y/K'* + AK') , (3) 

where K' = K(T)/C T - Since K'(T) is typically an in- 
creasing function of T, p(T) will generally decrease with 
increasing T. In our calculations, we have also assumed 
the simple van 't Hoff behavior 

K(T) = exp[-AG/k B T], (4) 

with a Gibbs free energy of formation 

AG(T) = Cl (T-T M )+c 3 (T-T M f, (5) 

choosing the values of c\, C3, and T M to be consistent 
with experiments on these DNA molecules [35] . The tem- 
perature T M can be interpreted as the melting tempera- 
ture of a single DNA duplex. 

Given p(T), l—p e fj is the probability that none of the 
duplexes forms a double strand. If the duplexes react in- 
dependently, this probability is simply 

l-p e// (T) = [l-p(T)] JV -/«, (6) 

where is the number of single strands on each 
monomer, and z is the number of nearest neighbors 
for the given lattice (z = 6 for a simple cubic lattice). 
Thus, p e ff(T) is also a decreasing function of T, but for 



iV,; > 1 will typically vary much more sharply with T 
than 1 -p(T). 

The criterion for the melting temperature Tm of the 
aggregate is easily written down for a fully occupied lat- 
tice. If we denote the melting probability by pm = 
p e ff(T M ), the melting temperature for a periodic lattice 
of monomers is simply 

l-p M = [l-p(T M )} Nd/z . (7) 

Eq. (7) implicitly determines T M in terms of Pm, N d} 
and z. As discussed in Ref. [18], if we assume that the 
specific links which are occupied at temperature T are 
time-independent, the melting takes place at pm = Pc, 
where p c is the bond percolation threshold for the lattice 
considered, at which an infinite connected path of double 
DNA strands first forms. For example, p c ~ 0.25 on a 
very large simple cubic lattice. 

Note that, according to eq. (7), p(Tm) decreases with 
increasing Nd and, hence, with increasing particle radius, 
since N d should be proportional to the surface area of a 
nanoparticle. To obtain specific values for Tm, we as- 
sume N d oc a 2 , set z = 6 and use the experimental result 
that N d = 160 when a = 8 nm [36]. Since p(T) de- 
creases monotonically with T, T M should thus increase 
monotonically with a, as reported in experiments [20,22]. 

In the present paper, we keep our structural model as 
simple as possible, so as to focus on the essential features 
of the aggregation process. In reality, there are many 
complicated issues which should be considered. For ex- 
ample, the binding of DNA to the gold nanoparticles is 
a statistical process [7], and thus there will always be a 
distribution in the number of DNA strands per particle. 
Moreover, because of the high local dielectric function 
arising from the large amount of DNA on the particles, 
not every DNA single strands which is bound on the 
particles can hybridize with the linker DNA [36]. Fur- 
thermore, it is unlikely that the gold/DNA composite 
resulting from the aggregation process will form a regu- 
lar crystal. All these features could, however, be included 
in an extension of our theory. 

III. OPTICAL PROPERTIES 

A. Discrete Dipole Approximation 

Given the distribution of clusters, and their geome- 
tries, we calculate their optical properties as follows. We 
first assume that, at any given time t, the extinction co- 
efficient C ext (t) of a given cluster can be calculated as if 
none of the other clusters were present. This amounts to 
neglecting corrections due to multiple scattering among 
the different clusters. We calculate C ex t for each cluster 
using the so-called discrete dipole approximation (DDA) , 
first proposed by Purcell and Pennypacker [37] . 
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As originally formulated, the DDA permits one to cal- 
culate the extinction coefficient of an irregularly shaped 
object of complex, frequency-dependent dielectric con- 
stant e(w), embedded in a homogeneous medium of 
real dielectric constant eh and subjected to an applied 
electromagnetic wave with electric field E cxp(k • r — 
ut) = Eo(r)e _mt . (In our notation, the physical field 
Ephj /S (r, t) is the real part of this complex quantity.) In 
the DDA, the object is replaced by a collection of N par 
identical point objects with polarizability a{ui) placed on 
a simple cubic lattice having a suitable lattice constant 
d. 

The relation between the polarizability a of these point 
objects and e(w) is discussed in Section IIIB. The induced 
dipole moment Pi of the i th polarizable point object is 
expressed as 



p l = qE Io 



(8) 



where E/ OCii , the local electric field at the position of the 



:th 



J loc,i 

point dipole, is 



E 



loc,i 



= E exp[ik • — iu)t] 



Pj- 



(9) 



Thus Eioc.i is the sum of the applied field at and all 
the scattered fields — Ajj - Pj emanating from the induced 
dipoles at Vj. In the DDA, the product A^ ■ pj can be 
expressed as [15,37,38] 



ikra —iwt 



-{k 2 r tJ x [( r<j x Pj ] + 



+ 



1 — ikr. 



r 2 L' ij 

ij 



fojPj -3r tf (r tf -Pj)]}. (10) 



Here = Vi — Tj and k = ui/c = 2n/ A, c being the speed 
of light in vacuum, and A the wavelength in vacuum. Eqs. 
(8)-(10) form a coupled set of 3N par equations, which 
can be solved for the N par dipole moments Pi using the 
complex-conjugate gradient method combined with fast 
Fourier transforms [38]. Given the p^'s, C ext for a given 
cluster is obtained from the relation 



Cext — 



4irk 
Wo 



-Im 



V ESexpHkTjO-pj 



(11) 



where the sum runs over the N par particles in the cluster. 



B. Dipole Polarizability 

In the DDA, there are many possible choices for the 
relation between the polarizability a(w) of the point ob- 
jects and e(oj) [39]. One possible choice is to use the 
Clausius-Mossotti equation: 



EM 



1 - (47r/3)na' 



(12) 



where n = l/d 3 is the number of point polarizable ob- 
jects per unit volume. This relationship between e(uj) 
and a is accurate for point polarizable objects on a cubic 
mesh, provided that the wavelength in the medium, A m 
is much larger than d [40]. However, if d/X m is not very 
small, this choice can violate the optical theorem. To pre- 
vent this violation, an extension of the Clausius-Mossotti 
equation to include a radiative reaction correction has 
been proposed [41]. For a simple cubic lattice, this cor- 
rection can be incorporated by using a lattice dispersion 
relation, which is appropriate for a periodic cluster [42]. 

One can also choose a by connecting it to the first 
scattering coefficient in the Mie theory, usually denoted 
a\. If we take the "point dipole" as a sphere of radius 
a, dielectric constant e, in a host medium of dielectric 
constant eh, then, as discussed in Ref. [43], the relation 
between a and a\ is 



Here 



ai = 



mil'i(mx)'ip[(x) — ipi(x)tp[(mx) 
rmpi(mx)^' 1 (x) — £i(x)ip[(mx) 



(13) 



(14) 



where ipi(x) — xj\(x) and £i = xh^\x), the complex 
number m = yje{ui)/eh, and x = 2ir^/ea/\. Here ji(x) 



(i) 



is the usual first order spherical Besscl function, and h\ 
is the first order spherical Hankel function. 



C. Application of the DDA to the Gold/DNA 
Composite System 

In the present work, each cluster consists of a number 
of DNA-linked individual gold nanoparticles. The ap- 
plication of DDA to this gold/DNA system clusters has 
been developed extensively for regular clusters [15,16]. In 
carrying out the calculations for irregular clusters, we do 
not explicitly include the optical properties of the DNA 
molecules, since these absorb primarily in the ultravio- 
let [17]. We use tabulated values for the complex index 
of refraction nb u ik{u) of bulk gold [44,45], then calcu- 
late C ext for each cluster using the DDA, using a finite- 
particle-size corrected dielectric function for the gold par- 
ticles," denoted e(u). 

We obtain e(w) by correcting the bulk dielectric func- 
tion ei,uik(w) = nl u[k (u)) to account for the additional 
damping mechanism induced by collision of the conduc- 
tion electrons with the particle surface. Specifically, we 
write 



e(u) = e bulk (uj) + 



ljj(uj + i/T) LU + i/T + i/Ta' 



(15) 
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where tbuik{u), lu p , and 1/t are the experimental bulk 
metal values for the dielectric function, plasma angular 
frequency, and relaxation rate. The quantity 



show the radius of gyration R g for such a cluster, plotted 
as a function of N m . R g is defined by the relation 



1/to = Av F /a eff 



(16) 



is the surface damping term, vp is the Fermi velocity, 
and a e fj is an effective particle radius, defined by set- 
ting 47rOgjy/3 = Vp ar t, where v part is the particle vol- 
ume. (Obviously, this expression is exact for spherical 
particles.) The constants uj p , t, and vf are taken from 
Ref. [46]. The constant A is a theory-dependent param- 
eter that includes details of the scattering process [47], 
and is expected to be of order unity. It also can depend 
on "chemical interface damping" (i. e. transfer of sur- 
face plasmon energy into excitation modes of the surface 
metal-adsorbate complex) [48] . To obtain the parameter 
A, we compared our calculations in the dilute limit to ex- 
perimental data for the extinction spectra of DNA-linked 
gold dispersed colloids [9], choosing the parameter A as 
the best fitted value, as further discussed below. 

Since we are treating disordered clusters, we believe 
that the best choice for a is eq. (13), and have used 
this equation in our calculations. For this choice, the 
DDA calculation has been compared to a more accu- 
rate method [49] for calculating the extinction coeffi- 
cient for a compact spherical aggregate of 89 30-nm gold 
nanospheres in an aqueous medium. The DDA result was 
found to be reasonably consistent with this more accurate 
method [16]. Also, for a 40-sphere DLA fractal cluster, 
Mackowski [50] has compared the exact total cross sec- 
tion, as obtained using a multiple scattering formalism, 
to that obtained from the dipole method in combination 
with cq. (13), using various choices for the sphere index 
of refraction. The calculated DDA results agree qual- 
itatively with the exact multiple-scattering calculation, 
except for a single discrepancy, which can be corrected 
by adding a frequency-independent constant to the ex- 
tinction coefficient. 

To improve the statistics, we average the calculated 
C ext for each cluster over possible orientations. We then 
sum the averaged extinction coefficients of all the indi- 
vidual clusters to obtain the total extinction coefficient of 
the suspension. This method should be adequate so long 
as the total volume fraction of clusters in the suspension 
is small (dilute regime). 



IV. NUMERICAL RESULTS 

We turn now to our numerical results, based on this ap- 
proach to modeling the structural and optical properties 
of gold/DNA composites. We begin with the structural 
properties. 

First, we show that our numerical algorithm does in- 
deed generate an RLCA cluster at p e ff = 1. In Fig. 1, we 



Rl 



1 



N„ 



Ei"- ! 



(17) 



i=l 



where is the position of the i th monomer and r = 
A" 1 J2 i rj is the cluster center of mass, df is then given 
by 



dlnN n 
hm — 

ii„^oo d In R„ 



(18) 



Fig. 1 shows that the log-log plot of N m against R g is 
indeed a straight line with a slope df = 2.1, consistent 
with expectations for RLCA, provided N m exceeds about 
200. 

Next, we consider how the radius of gyration R g , and 
hence the fractal dimension of the clusters, varies with 
Nmc- In Fig. 2, we plot R g as a function of the number 
of MC steps Nmc, for probability p e ff = 0.9. We con- 
sider the evolution of clusters having N m varying from 
100 to 600. In several cases, the final result of the evolu- 
tion is more than one cluster; in these cases, we calculate 
R g for the largest cluster. As shown in Fig. 2, as the re- 
structuring of cluster proceeds with increasing Nmc, the 
radius of gyration R g of the largest cluster of the system 
becomes smaller and eventually relaxes to a saturated 
value. 

To characterize these changes, we calculate the frac- 
tal dimension df at the beginning of the simulation, and 
after the system has fully relaxed. In Fig. 3, we show 
log- log plots of R g (N m ) for two different values of Nmc, 
using the data of Fig. 2. For small N M c {Nmc = 100), 
we measure df ~ 2.1, as expected for an RLCA cluster. 
But when we measure the value df for the saturated R g 
at large Nmc, we find df ~ 3, corresponding to com- 
pact, non-fractal clusters. Thus, the clusters are becom- 
ing more compact with increasing simulation time. 

In Fig. 4, we show some typical cluster morphologies 
corresponding to the procedure described above. In Fig. 
4(a)-(c), we show a cluster with N m = 1000, at p = 0.9, 
after 0, 7000, and 70000 MC steps. The gradual transi- 
tion from a fractal morphology to a more compact one 
is evident in the figure. (This behavior can also be ob- 
served at p = 0.95 although the saturation time is longer 
than at p = 0.90.) Also during the simulation, we ob- 
serve that most configurations have one large cluster and 
a few monomers as can be seen in Fig. 4 (b) and (c) . This 
indicates that at large p e f / the mechanism for cluster re- 
structuring is mainly the diffusion of monomers along the 
surface of the large cluster. 

From these numerical results we may infer some quali- 
tative conclusions about the sol-gel transition mentioned 
earlier. We denote by t\ the time needed to break one 
monomer from the surface of a large cluster. We esti- 
mate this time as t\ oc t /(1 — p e ff) Nav , where tq is the 



6 



time needed to break a single link between two particles, 
and N av is the average number of links on one monomer 
at the surface of a large cluster. Thus n diverges as 
Peff — * 1) but is finite for any p e ff < 1, allowing the 
largest cluster to relax to a non-fractal structure. Hence, 
we can expect that psa should be close to 1, if sufficient 
relaxation time is provided. 

Next, we turn to the calculated optical properties of 
these gold nanoparticle/DNA composites, based on this 
structural model. We begin by showing in Fig. 5 the 
computed extinction coefficients as a function of A for 
gold monomers with 6.5 nm radius, including the quan- 
tum size corrections embodied in eqs. (11), (13), (15), 
and (16). For comparison, we also show C ext (X) with 
no quantum size correction, and the experimental data 
in Storhoff, et al. [9]. The rather sharp extinction peak 
near A = 530 nm corresponds to the well-known surface 
plasmon peak, in which light is absorbed by an oscilla- 
tion of the electronic charge within the gold nanoparticlc. 
We find that the quantum-size-corrected extinction coef- 
ficients give the best fit to the experimental data if we 
choose A — 0.85. This value is of order unity, as ex- 
pected. 

We now discuss the optical properties of various clus- 
ters. In Fig. 6, we show the specific extinction coefficient 
Cext (A)/L 3 , as calculated using the DDA for clusters of 
gold monomers of radius 20 nm, arranged on a simple 
cubic lattice of edge varying from L = 1 to L = 7 lattice 
constants, and thus from N m = 1 to N m = 343. We 
choose the lattice constant £ — 48 nm. For L = 1, the 
extinction coefficient corresponds to a single monomer. 
As L increases, the peak first shifts towards the infrared, 
then broadens substantially. 

In Fig. 7, we show C ex t(X)/N m as a function of A for 
N m varying from 1 to 343, but with particles now form- 
ing a single RLCA cluster. Once again, the particles have 
radius 20 nm as in Fig. 6, and the RLCA cluster forms 
on a simple cubic lattice of the same lattice constant 
I = 48 nm. The cluster is generated using the RLCA 
algorithm described at the beginning of Sec. IIB. As in 
Fig. 6, the extinction peak is red-shifted and broadened 
as N increases, but both the shift and the broadening are 
substantially smaller than in Fig. 6. 

As a final comparison, we show in Fig. 8 the specific 
extinction coefficients C ex t(X)/N m for the fully relaxed 
configurations generated by the algorithm of Section IIB. 
We set p e f f = 0.9 and consider N m varying from 1 to 400, 
as indicated in the legend. The behavior of C ex t in this 
case resembles Fig. 6 more than Fig. 7; in particular, the 
surface plasmon peak is more broadened and red-shifted 
than that of Fig. 7. This result is not surprising, since 
the clusters of Fig. 8, like those of Fig. 6 and unlike those 
of Fig. 7, are compact and non-fractal, though they in- 
corporate some disorder. 

Figs. 6 and 8 show that compact, non-fractal clusters 
(whether regular or irregular) have similar optical prop- 



erties. Now, these clusters do have various structural dif- 
ferences; for example, the clusters of Fig. 8 have rougher 
surfaces than those of Fig. 6. We conclude that these 
structural differences do not significantly affect the clus- 
ter extinction coefficients. Of the three extinction co- 
efficients shown in Figs. 6-8, that of Fig. 7 agrees best 
with experiment; those of Figs. 6 and 8 show too large a 
broadening and red shift of the surface plasmon peak. 

In Fig. 9, we show the calculated normalized extinc- 
tion coefficient C ex t(X)/N m at fixed wavelength A = 520 
nm as a function of p e f / . The three continuous curves 
represent the results of the model described in our previ- 
ous paper [18]; In this model, the cluster is constructed 
starting from a simple cubic lattice of linear dimension 
Li, fully occupied by monomers of radius a. Bonds are 
then removed with probability 1 — p e ff ■ At a given value 
of Peff, there are one or more clusters, depending on the 
relation of p e f / to the bond percolation threshold p c for 
the lattice and on L. The open circles represent the re- 
sult of the model of this paper - that is, the final cluster 
is the result of a series of MC steps, each consisting of 
an RLCA aggregation, followed by a thermally induced 
breakup of the cluster. The calculated C ext {X)/N m are 
based on the long-time and thus relaxed configurations 
generated by this procedure. 

As is evident from Fig. 9, there is a sharp drop in 
Cext (A) /N m at a characteristic value of p e f / in both mod- 
els, which we identify with the melting transition for that 
model. For the bond percolation model of Rcf. [18], the 
melting value pu is close to 0.25, the bond percolation 
threshold on a simple cubic lattice in d = 3. For the 
model of the present paper, the melting point pm occurs 
at a considerably higher value of p e ff- In the present 
model, unlike the previous model of Ref. [18], the final 
cluster results from many steps of link removal and ref- 
ormation. Thus the infinite cluster which is formed from 
the percolation model, can not survive at the bond per- 
colation threshold because it contains many weak links 
which will break apart after several link removal steps. 
Thus, there should be an increase in the melting value of 
Pm, as observed in our simulations. 

In Fig. 4 (d) , we illustrate a typical configuration gen- 
erated at Peff = 0.6 by the present model. Here we can 
clearly see that there is no large cluster, even though the 
Pef f is much higher than the bond percolation threshold 

Pe- 

Even though the value of pu at melting is larger for 
the present model than in our previous model in Ref. 
[18], the sharpness of the melting transition (as seen in 
the calculated optical properties) is similar to that found 
previously. The reason is that the sharp transition re- 
sults mainly from the individual link properties [cf. eq. 
(3)], not the behavior of the cluster. Thus, the present 
model preserves the sharp melting transition reported in 
experiments [20,22]. 

In Fig. 10, we replot the results of Fig. 9, but with 
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p e f f translated into a temperature using the prescription 
of eqs. (2) and (3). The higher melting value of p e ff ob- 
tained from the present model now translates into a lower 
melting temperature Tm, in comparison to the percola- 
tion model in Ref. [18]. 

A striking feature of Figs. 9 and 10 is the "rebound 
effect," seen in the cross and square of Fig. 9, and in the 
cross in Fig. 10. The square at p e ff = 1 in Fig. 9 rep- 
resents the calculated C ext (X)/N m resulting from pure 
RLCA. The cross at p eff = 0.9 in Fig. 9 or in Fig. 10 
represents a configuration obtained after 7000 MC steps, 
which is not long enough to produce the asymptotic com- 
pact cluster for this process at this value of p e ff- In 
both cases, C ext {X)/N m is larger than that of the fully 
relaxed configuration for this p e ff- This behavior can 
be understood quite simply. Since these points are not 
fully relaxed, they correspond to clusters which are more 
fractal than the fully relaxed clusters. Being more frac- 
tal, they show some characteristics of the melted clusters 
at this wavelength; hence, they have a somewhat higher 
C ex t{X)/N m than the fully relaxed samples. In fact, this 
rebound effect has been observed in experiments which 
are carried out on gold/DNA aggregates [4,9,17,20,22]. 
However, when gold/DNA composites are formed in the 
presence of a DNA-coated flat surface, this effect was not 
observed [20]. 

The observation of the rebound effect in some experi- 
ments but not in others needs to be understood, as does 
our result that among the three extinction coefficients 
shown in Figs. 6-8, that of Fig. 7, which represents the 
fractal aggregate, agrees best with experiment. Both re- 
sults suggests that fractal-like gel state persists up to sur- 
prisingly high temperatures, such as room temperature. 
On the other hand, our numerical study of our structural 
model suggests that psa is either 1 or at least close to 
1, given a sufficiently long relaxation time, which implies 
that a gel state should exist only a very low tempera- 
tures at such long times. We will discuss both of these 
questions in the next section. 

V. DISCUSSION AND CONCLUSIONS 

In this paper, we have described a structural model 
for the formation of gold/DNA nanocomposites, going 
well beyond our earlier, purely percolation model. The 
model includes several key features expected to play a 
role in these aggregates. For example, we include both 
the formation of clusters via reaction-limited cluster- 
cluster aggregation, and the thermally induced dehy- 
bridization which breaks up the links between the gold 
monomers. We also calculate the aggregate extinction 
coefficient C ext {\) for some of the model structures. 

Our structural model is characterized by two types 
of transitions. The first is the melting transition men- 
tioned above. Above pm, in the limit of a large number 



of monomers, the system is characterized by at least one 
large cluster; below pm, this cluster breaks up into two or 
more smaller clusters. In our previous percolation model, 
Pm coincides with the percolation threshold, but in the 
present model, it occurs at a higher value of p e ff- It is 
observable as a rather sharp transition in C ext (X) at a 
characteristic wavelength of about 520 nm. The other 
transition is the sol-gel transition. For p > psc, the ag- 
gregate has a fractal gel-like structure, while for p < psa, 
it is a non-fractal sol. We have presented some numerical 
evidence showing that this sol-gel transition will occur at 
p e ff ~ 1, if the sufficient relaxation time is provided. 

However, our calculated C ex t(X) agree best with ex- 
periment if we assume that the composites grow into a 
fractal gel-like structure, rather than the compact sol 
structure. Specifically, the surface plasmon peak seen 
in C ex t{\) for isolated gold monomers is broadened and 
red-shifted by about the right amount, in comparison to 
experiment, for the fractal clusters, but by far too much 
for the compact, non-fractal clusters. 

Thus, it appears that the full restructuring of our 
model, which would eventually give rise to a non-fractal, 
does not go to completion in the experimental compos- 
ites. There are various possible reasons for this behavior. 
Most likely, the time required for a full relaxation is sim- 
ply too long. If the composite is well below its melting 
temperature, only a small fraction of links are broken and 
the relaxation is likely to be very slow. As an illustra- 
tion, we found that, for p e ff = 0.9, even Nmc = 7000 is 
not long enough to produce the relaxation necessary to 
generate a non-fractal cluster. 

Another suggestive piece of experimental evidence 
is the rebound effect. According to our calculations, 
CextW should show a characteristic increase at A = 520 
nm, above its fully relaxed value, when the cluster is 
relaxed for only Nmc — 7000 steps. With this in- 
crease, C ex t(X) is slightly larger than the fully relaxed 
value expected at this A for a non-fractal cluster. Ex- 
periments do indeed show indications of this rebound ef- 
fect [4,9,17,20,22]. 

In summary, we have presented a model for the struc- 
tural evolution of gold/DNA composites. This model 
leads to a wide range of structures which depend on 
the characteristic parameters of the model: the num- 
ber of monomers, the simulation time, and the temper- 
ature. The optical properties of these composites cal- 
culated from the resulting structure are consistent with 
experiment. In particular, we obtain the observed sharp 
melting transition with a characteristic dependence on 
monomer radius, a characteristic shift and broadening of 
the surface plasmon peak in the extinction, of about the 
experimental magnitude, and a characteristic rebound in 
Ce£ct(A), which indicates that the restructuring of the ex- 
perimental composites has not been able to go to com- 
pletion. In view of these successes, we plan to extend the 
present work to study the sol-gel transition in gold/DNA 
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nanocomposites. The results of this investigation will be 
described in a future publication [51]. 
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FIG. 1. Log-log plot of the number of monomers N m for an 
RLCA cluster at p e ff = 1, plotted as a function of the radius 
of gyration R g . Squares: present calculations; dashed line: 
least squares fit. The slope of the dashed line is 2.1 = df, 
consistent with expectations for RLCA clusters in d — 3. 
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FIG 2. Variation of R g with number of Monte Carlo steps 
Nmc, as obtained using the algorithm described in the text 
for probability p e ff = 0.9, and clusters with a number of 
monomers N m varying from 100 to 600, as indicated. 
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FIG 3. Log-log plot of the number of monomers N m for 
the clusters generated by the algorithm described in the text, 
at p e ff — 0.9, as a function of the radius of gyration R g . The 
open circles represent the data of Fig. 2 for N m as a func- 
tion of R g at Nmc = 100. The slope of dashed line is 2.1, 
the expected value for RLCA. The open squares represent the 
relaxed values of R g for large Nmc, also taken from Fig. 2. 
They are well fitted by the solid line, which has a slope of 3.0, 
corresponding to a compact, non-fractal cluster. 




FIG. 4. (a) Initial form of the cluster, after the first reac- 
tion-limited aggregation (Nmc = 0), for N m = 1000, p = 0.9. 
The clusters are divided from an RLCA fractal, with frac- 
tal dimension df ~ 2.1. (b). Same cluster, but after 7000 
Monte Carlo steps, (c) Same as (b), but after 70000 Monte 
Carlo steps (nearly saturated), (d). Clusters with p = 0.6, 
N m = 1000, and N M c = 6000. 
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FIG. 5. Extinction coefficient C ext {\) per unit volume of 
gold for a dilute suspension of gold nanoparticles in water, 
plotted for a particle radius a — 6.5 nm. Crosses: experimen- 
tal data of Ref. [9]. Full curve: calculated extinction coeffi- 
cient including quantum size corrections, as computed from 
eqs. (15) and (16) with A — 0.85. Dashed line: calculated 
extinction coefficient without quantum size corrections. 
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FIG. 7. C ext {\) I 'N m plotted as a function of A for various 
number of particles N m from 1 to 343. In all cases, we assume 
the cluster forms a RLCA cluster with particles of radius 20 
nm, and is generated on a lattice with lattice constant I — 48 
nm. Ceit (A) is calculated using the DDA as in Fig. 5. 
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FIG. 6. C elt (A)/L 3 plotted as a function of A for Lx Lx L 
clusters with edge varying from L = 1 to L = 7 lattice con- 
stants. In all cases, we assume the cluster forms a cubic com- 
pact cluster with particles of radius 20 nm, lattice constant 
I — 48 nm; so the number N m of monomers varies from 1 
to 343. We use the DDA including quantum size corrections 
from eqs. (15) and (16) and polarizability obtained from eqs. 
(10), (11), (13), and (14) 
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FIG. 8. Ce^t (A) /N m plotted as a function of A for a number 
N m of monomers varying from 1 to 400. In all cases, we use 
relaxed configurations from the algorithm in the text, whose 
radii of gyrations are shown in Fig. 2. (Examples of these 
configurations are shown in Fig. 4.) Once again, we assume 
a monomer radius is 20 nm, a lattice constant of 48 nm, and 
use the same algorithm as in Figs. 5 and 6 to calculate the 
extinction coefficient. 
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FIG. 9. Normalized C ext (\) at a fixed wavelength of 
A = 520 nm, as calculated in the DDA, for clusters of var- 
ious geometries and different numbers iV m of monomers of 
radius 20 nm, plotted as a function of p e ff- Full curves: clus- 
ters generated by percolation model of Ref. [18]; p e ff repre- 
sents the fraction of links which are present in this model. 
The total number of monomers is L 3 . Open circles: C ex t at 
A = 520 nm for the model of the present paper, with N m — 
1000, with fully relaxed long-time configurations. The square 
at p e ff = 1.0 represents the calculated C ex t at A = 520 nm 
assuming RLCA clusters. Cross at p e // = 0.9 represent the 
calculated values of C ex t at A = 520 nm for the unsaturated 
configuration obtained after Nmc = 7000 steps; the cross at 
p e ff = 1 is the value of C ex t obtained after a single RLCA 
aggregation. 
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FIG. 10. Same as Fig. 8, but with p e ff translated into a 
temperature, using eqs. (2) and (3). The cross corresponds to 
Peff = 0.9 after Nmc = 7000. 
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